A random matrix-like model for the Polyakov loop and center symmetry 
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We formulate a random matrix-like model for the Polyakov loop in SU{N) Yang-Mills theories. 
It describes a simplified dynamics in terms of eigenvalue differences. The deconfinement phase 
transition encoded in center symmetry breaking is reproduced including its nature being first order 
for SU{3) and second order for SU{2). Analytical arguments about the phases are presented and a 
comparison to other approaches is made. 
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Introduction — Evidence for the quark-gluon plasma 
state from heavy ion experiments has triggered plenty of 
theoretical and numerical work on the phase diagram of 
Quantum Chromodynamics (QCD). Deconfinement and 
restoration of chiral symmetry are the principal phenom- 
ena to be described at the phase transition. As intrinsi- 
cally nonperturbative effects they remain difficult to de- 
rive from first principles. In perturbation theory, for in- 
stance, quarks arc not confined. 

Random matrix theory - in contrast to ab initio- 
methods like e.g. lattice gauge theory - governs universal 
spectral properties. In QCD this usually concerns the 
Dirac operator. Through the density of Dirac eigenval- 
ues at zero eigenvalue, random matrix theories contain a 
chiral condensate quite naturally. Furthermore, temper- 
ature can be built into the formalism via the Matsuba- 
ra frequencies (see [1] and references therein). Although 
such random matrix models are schematic, they correctly 
incorporate the chiral restoration at high temperatures. 
In this way random matrix models help to understand 
one of the principal phenomena of QCD, in a very simple 
system with almost all QCD dynamics removed. 

In this Letter we ask the question whether in the same 
spirit one can investigate the finite temperature decon- 
finement transition, too. The Polyakov loop is related to 
the free energy of (infinitely heavy) quarks and hence the 
deconfinement order parameter at finite temperature. It 
can be studied in the quenched approximation aka Yang- 
Mills theories. In the latter case, the Polyakov loop is a 
strict order parameter as it respects/breaks center sym- 
metry at low/high temperatures, respectively (see be- 
low). 

We shall introduce a random matrix-inspired model 
for the Polyakov loop eigenvalues and their effective 
potential depending on temperature. As argued this 
is achieved already in pure Yang-Mills theories, with 
general gauge group SU{N). Our model reveals center 
phase transitions with the correct order for SU{2) 
and SU{3). Thus also center symmetry breaking and 
hence deconfinement can be understood by virtue of a 
simplified model based on random matrices. 

The Polyakov loop and center symmetry — The 
Polyakov loop as an element of the gauge group SU{N) 



can be diagonalized to 



V = exp(27ridiag(^i,^2,---,A'A')) 



(1) 



where the parameters fj,a are dimensionless numbers and 
sum to an integer. Of course, the /i's are defined only 
up to permutations and adding integers, which will be 
discussed in more detail below. 
The average over the trace 



1 1 ^ 



(2) 



is an order parameter for (de) confinement via center sym- 
metry [2] (see [3] for a recent review): The Polyakov 
loop is traceless, (L) = 0, below the critical temperature, 
which is achieved by equidistant eigenvalues 



fii = Afi, /i2 = A^ 



1 



, ^AT = A^ 



iV- 



N 



(3) 



where an additional shift A/i = 1/2N is needed for even 
N to ensure detT' = 1. 

At high temperatures the absolute value of the 
Polyakov becomes maximal, (|L|) — >■ 1, which amounts 
to degenerate eigenvalues 



fJ-i, M2, M3, 



,fj.N^ki-- 



(4) 



Here /x is for real Polyakov loop L — > 1, in the quenched 
case at high temperature ^ also approaches (with same 
probability) multiples of 1/iV with P being a nontrivial 
center element of SU{N) and L being a nontrivial iVth 
root of unity. 

This reflects center symmetry, which is the invariance 
under transformations 



L -> exp(27riA:/A^)L , 
with integer k. 



fJ.a^ kl-a + 



N 



(5) 



"Deriving" the model — As is natural for random ma- 
trix models, we will consider constant fields. The tempo- 
ral gauge field ^o needs to reproduce the Polyakov loop, 

^0 = T-jlogP = 27rTdiag(^i,^2,---,A'Jv)- (6) 



This choice amounts to the Polyakov gauge whose 
Faddeev-Popov determinant (Jacobian) is known to be 
the reduced Haar measure [4] 



Y[ sin^ (t^CMo - M&)) 



(7) 



a>b 



This kinetic term gives sets of eigenvalues fj,a different 
weights according to the associated volume of gauge or- 
bits. 

In our model the spatial gauge fields Ai will take the 
role of fluctuations and hence are replaced by hermitian 
matrices Xi. 

The action proportional to F^^ = X^u i^ t'^C^I^A" "^i^])^ 
is treated as follows: The magnetic part 

— J2i tr[^n Aj]"^ is replaced by the quadratic J2i tr^j^j 
Gaussian weights for the random matrices Xi . Assuming 
universality, the precise form of the weight for Xi should 
not influence the results much. 

The electric part — J2i ^A-^o, AiY is quadratic in both 
Ai and A^. Upon replacing the gauge fields by the ma- 
trices Xi and the rhs. of Eqn. (pi), respectively, this part 
consists of 2 X 2-blocks 



-tr 



Ma 



Mb 



xf* 



xf 
xf 



ab\2 



2{^la-^,,Y\xf\ 



(8) 
to be muhiplicd by (27rT)2. 

Here the ambiguity in the /i's leads to a subtlety. 
Consider momentarily ordered /i's according to (e.g. [5]) 
Ml ^ M2 ^ ■ ■ ■ < Mw < Ml + 1 — Ml where the equiva- 
lence = stands for the same eigenvalue exp(27ri/ia) in the 
Polyakov loop. Then for two indices a > b: 



^J■b < ^J■a < fj-b + 'i- = ^J■b ■ 



(9) 



The block considered so far, Eqn. ([8]), takes into ac- 
count the commutator of the spatial gauge field with 
Aq ^ diag(Ma, fJ-b) resulting in a term with the difference 
^J'a — Mfc- However, Aq ^ diag(Ma,M& + 1) is equally valid 
and so is the difference Mb + 1 — Mo = 1 — (Mo — fJ-b)- The 
above inequalities limit these differences to the interval 
[0, 1] and it depends on the configuration which one is 
shorter (yielding a smaller action). 

For symmetry reasons we therefore include also the 
other difference [12], through a term 



fab\2 



(i-(Ma-Mb)) vr'i 



(10) 



with new fluctuation matrices X^ subject to the same 
quadratic action X^i^rX^^. 

Both X and X' are N x A^-matrices for gauge group 
SU{N) (unlike in conventional random matrix models, 
where the limit of infinite size is performed). 

As our model depends only on the differences of m's 
center symmetry is manifest. We will abbreviate the dif- 
ference of two m's by V. The ordering of the eigenvalues 
can be taken into account by the definition 



in terms of which we can finally write down our model: 
Z({/i}) = Y[ sin^ {niyab) (12) 



a>b 



Yl fdX,dXlexp(- (27rC)2tr(A,2 + Af ) 

exp(- ^2(2^r)2[^„\|Af p + (1 - i^abflX',^']'] 



a>b 

thereby introducing the parameter C, a 'coupling con- 
stant' of the model. 

The effective potential for the Polyakov loop eigenval- 
ues is related to Z hy Z — exp{—V/T), hence up to 
z^-independent constants 



^^ = ^{ -log sin^ (...,) (13) 

a>b 

with Ng the number of spatial dimensions {i — 1, . . . , Ng 
in ( 12 )), in due course put to 3 unless specified otherwise. 



Results and interpretation — In the dynamics of 
Polyakov loop eigenvalues the Haar measure term sup- 
presses degenerate eigenvalues, for which i^ab = 0, and 
hence tends to confine, as is known from strong coupling 
expansion on the lattice, see [5] and references therein. 

The I'^j-term in the second line of Eqn. (13) prefers 
nearby eigenvalues, which however renders the (1 — i^ab)^- 
term large (for fJ-a — l^b ^ £; the terms act vice versa for 
tJ-a = Mb ~ f ) ■ Which ij, configurations eventually yield the 
minimal potential V then depends on the accompanying 
factors and hence on temperature measured in units of C. 




^ab = [tJ-a - i^b) mod 1 G [0, 1) 



(11) 



FIG. 1: Effective potential V/T ove r the traced Polyakov l oop 
L for gauge group SU{2), Eqn. ( 14 1 (obtained from Eqns. | |13[ ) 
and ([2|), at low temperature T/C =1.0 (upper panel) and 
high temperature T/C — 6.0 (lower panel), respectively. 
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FIG. 2: The inverse effective potential —V/T over the traced 
Polyakov loop L for gauge group SU(3), Eqns. ( |13[ ) and (pi, 
at low temperature T/C = 3.0 (upper panel) and high tem- 
perature T/C = 6.5 (lower panel), respectively. 



In Fig. [T] we show for two temperatures the effective 
potential in SU{2) over the traced Polyakov loop L as 
given by 



2T 



= -log(l-L2 



(14) 






^arccosLN^ 2\ fC^ 

Vr2 



/ arccos L^2 



The minimum of the effective potential moves from i = 
at low temperatures to \L\ -^ 1 at higher ones as it 
should. The invariance under L — > —L reflects the man- 
ifest center symmetry of the effective potential. 

The corresponding plots for gauge group SU{3) are 
shown in Fig. [2] Here V is plotted over the com- 
plex L, whose values are bounded by an approximate 
triangle. Again, L = is preferred at low tempera- 
tures, whereas at high temperatures the minima are at 
L — > {1, ±exp(27ri/3)} as expected. 

We note that for SU{2) and SU{3) there is a 
one-to-one mapping of Polyakov loop traces L (real 
resp. complex) to configurations of eigenvalues fia (one 
resp. two independent ones). A fixed L in SU{N) 
with iV > 4 is generated by different /i^ 's and hence has 
different potentials V (the minimal one should be taken) . 



Our model correctly predicts the confining, center- 
respecting Polyakov loop L = at low temperatures and 
a breaking of center symmetry by the 'ground state' with 
\L\ -^ 1 in the high temperature deconfining phase (with 
L — turned into a local maximum). 

What is more, also the order of the transition in our 
model is the well-known one for Yang-Mills theories. As 
can be read off from Fig. |3) the phase transition is second 
order for SU{2) since the finite i-minimum develops in 
a continuous manner out of the L = 0-minimum. In 
contrast, the transition is first order for SU{3), since 
there is a phase of coexisting local minima, above which 
the finite L-minimum becomes the global one. 

The low and high temperature behavior of the effective 
potential can be understood through approximations of 
the second line of Eqn. (13 1. At low T the argument of 



the log is expanded around the corresponding power of 
C'^/T'^, the first j/-dependent term is proportional to the 



Ea 



>b 



3 ^ i^ab + (1 - '^abf) 



(15) 



This is maximal at Vab — and thus suppresses degener- 
ate eigenvalues like the Haar measure term. Because of 
the T^-prefactor in ( 15 1 the Haar measure is actually the 



leading term at low temperatures resulting in confining 
Polyakov loops. 

At high T the leading log-term is the sum '^^j^yi, over 



3\0g{{lyab)\l-'^abf)^ 



(16) 



As the argument is positive it is minimized by vanishing 
Vab, hence this term prefers deconfining degenerate eigen- 
values. Whether a deconfining Polyakov loop is taken on 
eventually depends on the dominance of this contribution 
over the Haar measure — logsin^(7rJ/Qb). As is clear from 
the figures [T] and [2] this is indeed the case for the physical 
situation of three spatial dimensions. 

For arbitrary number Ns of spatial dimensions the sec- 
ond z^afc-derivative at intermediate Vab = 1/2 of the term 
( 16 1 and the Haar measure is — 16A^s vs. 27r^, respectively. 



This suggests a transition to deconfinement at high tem- 
peratures for Ns > 2. This argument is exact for SU{2) 
~ where there is just one independent eigenvalue differ- 
ence 1/21 ~ but is supposed to hold also for higher gauge 
groups. 

Interestingly, the effective potential at high tempera- 
tures, Eqn. (16), has as argument the one-loop result [Ij. 
It follows that both are minimized by the same /i's. How- 
ever, our V is the log of the one-loop expression. There 
seems to be no reason why the random matrix model 
built on constant gauge fields could exactly reproduce 
perturbation theory or effective Polyakov loop models [7] 
which are based on varying gauge fields (and typically 
invoke character expansions). 




Summary and outlook — We have formulated a 
randoni-inatrix like model for the Polyakov loop eigenval- 
ues, Eqn. (13 1 (using continuum language). It contains a 



FIG. 3: Effective potential V/T (with V/T\l=o subtracted) 
as a function of L zoomed in for three temperatures near 
the phase transition. Upper panel: SU{2) with T/C = 
3.62, 3.637, 3.655, a second order transition. Lower panel; 
real Polyakov loops in SU{3) (parametrizing ni = —fig £ 
[0,1/2], ^2 = 0, L = [1 + 2cos(27rpi)]/3) with T/C = 
4.53, 4.5635, 4.588, a first order transition. 



Phenomenological approaches like Ginzburg-Landau 
for second order phase transitions or [5] for first order 
ones use low powers of the order parameter near the crit- 
ical temperature. We make the common ansatz 



"^ILP 



64, 



-(L^ + L^) + -^|Lr (17) 



^phen(-^) 

T ^ 2 '"' 6 '" ■ " ' 4 

where the cubic term is utilized only for SU{3). 

The coefficients 62,3,4 can be obtained in SU{2) by 
Taylor expanding Eqn. ( 14 ) around L = to fourth or- 



der. In SU{3) we compare the (^i,^2)-Taylor series of 
F(^i,^2,-Mi-M2) with VphcniL{^J'l,^J■2,-^H-^J■2)) near 
the confining value (/ii,^2) — (1/3,0). 

In our model, 62,3,4 come out as rational functions of 
T/C [not displayed]. In SU{2) the quadratic coefficient 
62 turns positive at a critical T/C = 3.637 (the solution 
of a certain quartic equation) signalling a second order 
phase transition, cf. Fig. [s] upper panel. In SU{3) this 
coefficient turns negative at T/C = 4.718. The finite L- 
minimum develops at a slightly lower critical temperature 
T/C ^ 4.56, cf. Fig. flower panel. 

Around these temperatures the higher coefficients 63 4 
are positive which is expected and e.g. agrees qualita- 
tively with the fits given in [3]. While the polynomial 
approximation to V describes the transition region fairly 
well, it breaks down at the latest at T/C > 4.88 resp. 
T/C > 6.51 where the quartic coefficients for SU{2) resp. 
SU{3) turn negative and higher powers of L become im- 
portant. 



deconfinement transition of correct order for SU{2) and 
SU{3). Through the (symmetrized) interaction with spa- 
tial gauge fields the Polyakov loop eigenvalues change 
from equidistant to nearly degenerate thereby breaking 
center symmetry. Hence this simple model captures the 
relevant finite temperature mechanism with only one pa- 
rameter C setting the temperature scale. 

The underlying space only enters through its dimen- 
sionality (the number of spatial gauge fields) hence the 
model neglects e.g. Polyakov loop clusters and topologi- 
cal objects that might drive the transition. 

In the quenched case the model could be used to inves- 
tigate the phase transition in the large TV limit (for the 
Wilson loop eigenvalues universal properties have been 
analyzed in this limit [10]) and in exotic gauge groups. 
In SU{3) there are also lattice data for the Polyakov loop 
in higher representations jll j to compare to. 

Furthermore, it should be straightforward to include 
quarks into the model, along the lines of the Dirac 
operator in chiral random matrix theory. We conjecture 
a weakening of the phase transition. In turn, our 
approach makes it possible to incorporate the Polyakov 
loop and hence center symmetry in random matrix 
models. Obviously an extension of the model to nonzero 
chemical potential would be interesting, too. 

The author likes to thank Jacques Bloch, Kenji 
Fukushima, Christof Gattringer, Tilo Wettig, Andreas 
Wipf and in particular Jac Verbaarschot for many help- 
ful discussions. This work has been supported by DFG 
(Contract No. BR 2872/4-2). 



[1] J. J. M. Verbaarschot and T. Wettig, Ann. Rev. Nucl. 
Part. Sci. 50, 343 (2000), hep-ph/0003017 

[2] A. M. Polyakov, Phys. Lett. B72, 477 (1978); 
L. Susskind, Phys. Rev. D20, 2610 (1979). 

[3] R. D. Pisarski, Nucl. Phys. Proc. Suppl. 195, 157 (2009). 

[4] N. Weiss, Phys. Rev. D24, 475 (1981); D. J. Gross, 
R. D. Pisarski, and L. G. Yaffe, Rev. Mod. Phys. 53, 43 
(1981)2; N. Weiss, Phys. Rev. D25, 2667 (1982). 

[5] T. C. Kraan and P. van Baal, Phys. Lett. B435, 389 
(1998), hep-th/9806034, D. Diakonov and V. Petrov, 
Phys. Rev. D76, 056001 (2007), 0704.3181. 

[6] K. Fukushima, Prog. Theor. Phys. Suppl. 153, 204 
(2004), hcp-ph/0312057, 

[7] M. Engelhardt and H. Reinhardt, Phys. Lett. B430, 
161 (1998), hep-th/9709115; A. Dumitru, Y. Hatta, 
J. Lenaghan, K. Orginos, and R. D. Pisarski, Phys. 
Rev. D70, 034511 (2004), hep-th/0311223, D. Diakonov 
and M. Oswald, Phys. Rev. D70, 105016 (2004), 
hep-ph/0403108; A. Wipf, T. Kastner, C. Wozar, and 
T. Heinzl, SIGMA 3, 006 (2007), hep-lat/0610043, 

[8] A. Straley and M. Fisher, J. Phys A6, 1310 (1973). 



[9] C. Ratti, M. A. Thaler, and W. Weise, Phys. Rev. D73, 

014019 (2006), hep-ph/0506234, [11] 

[10] D.J. Gross and E. Witten, Phys. Rev. D21, 446 

(1980); B. Durhuus and P. Olesen, Nucl. Phys. B184, [12] 

461 (1981); A. M. Brzoska, F. Lenz, J. W. Negele, 

and M. Thies, Phys. Rev. D71, 034008 (2005), hepj 



th/0412003 , R. Narayanan and H. Neuberger, JHEP 12, 



066 (2007), 0711.4551. 

S. Gupta, K. Hiibner, and O. Kaczmarek, Phys. Rev. 

D77, 034503 (2008), 0711.2251. 

One could also include higher differences such as /^a + 2 — 

Hb, but we want to keep the model as simple as possible. 



